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SENSITIVITY  ANALYSIS  HITS  A  SIMPLE  CHARGING  CODE 

by 

A.  J.  Sims 
C.  L.  Urenn 


SUMMARY 


Computer  simulation  techniques  using  codes  such  as  NASCAP  have  proved  very 
successful  in  tackling  the  spacecraft  charging  problem.  Although  a  full  dynamic 
three -dimensional  cods  is  required  to  model  a  spacecraft  in  orbit,  simpler  codes 
(eg  MATCHG)  can  do  much  to  clarify  the  interaction  bepveen  the  environment  and  a 
surface.  It  is  difficult  to  obtain  complete  data  on  either  the  incident  particle 
fluxes  or  the  material  properties  vhich  determine  the  various  current  components. 
The  latter  are  dependent  upon  particle  energy  which  is  in  turn  a  function  of 
surface  potential;  therefore  the  resulting  equilibrium  and  tha  charging  time  to 
attain  this  potential  are  controlled  by  a  long  list  of  input  parameters.  In  order 
to  assets  the  sansitivicy  of  a  solution  to  each  of  these  inputs,  a  simple  coda 
(EQUIPOT)  it  used  to  calculate  the  net  current  as  surface  potential  is  stepped 
from  zero  to  that  required  for  equilibrium.  jJnat  ■  The  required  accuracy  for 

all  the  parameters  it  chon  established  by  successively  introducing  a  delta  of  1 0X 
on  each.  This  sensitivity  analysis  is  performed  or.  the  properties  of  a  shaded 
kapeon  patch  on  a  spherical,  conducting,  sunlit  spacecraft  subject  to  an  ambient 
plasma  environment  vhich  will  produce  large,  negative  potentials..  Two  main  results 
are  obtained:  firstly,  identification  of  the  important  role  played)  by  the  high 
incident  energy  part  of  the  secondary  electron  yield  curvs  and  secondly  the 
ability  of  a  low  density,  cold  plasma  component  to  limit  negative  charging. 
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1  INTRODUCTION 

Much  effort  ha*  baan  expended  in  th*  development  of  software  u*ad  to  nodal 
the  charging  of  spacecraft  by  nagnetoapharic  plasmas  .  Fully  dynamic,  three- 
dioansional  coda*  ara  a**anr.ial  for  dacarmining  tha  complex  pocantial  contours 
and  diffarancial  pocantial*  which  davalop  around  a  ipacacraft  compoeed  of  many 
diffarcne  aurfaca  matarial*  in  a  cotnpl  icatad  gaomatry.  In  order  for  such  codes 
to  giva  reliable  results,  it  is  aasential  that  the  basic  interaction  between  the 
plasma  environment  and  a  given  surface  matarial  is  wall  understood. 

By  considering  incident  ion  and  alactron  currants  to  a  surface  (a  function 
of  thu  plasma  environment  and  surface  pocantial)  and  tha  secondary  currants 
emitted  from  a  surface  (a  function  of  the  incident  currants  and  surface  material 
properties)  it  is  possible  to  construct  a  simple  numerical  simulation  to 
datarmina  tha  equilibrium  surface  potential  at  which  the  net  current  to  the 
surface  is  saro.  Such  a  coda  may  ba  used  to  asses*  the  sensitivity  of  the 
solution  to  small  changes  in  the  input  parameters  which  define  the  plusma 
environment  and  tha  surface  matarial  properties.  The  results  of  such  a  sensi¬ 
tivity  analysis  provide  an  estimate  of  tha  error  bars  on  the  results  of  more 
cosg>lex  charging  codes  and  also  indicate  where  more  accurate  material  data  or 
plasma  measurements  must  be  made. 

Tha  EQUIPOT  code  has  been  developed  to  compute  the  equilibrium  potential 
of  a  Surface  in  a  given  plasma  environment,  and  therefore  provides  an  ideal 
framework  for  carrying  out  a  detailed  sensitivity  analysis.  This  Memorandum 
gives  a  brief  overview  of  tha  code  and  presents  th*  results  of  a  s*  sitivicy 
analysis  performed  using  the  code.  The  scope  of  the  analysis  described  here  is 
limited  to  the  study  of  a  shadowed  or  partially  illuminated  kapton  patch  on  a 
sunlit,  conducting  aluminium  spacecraft  structure.  Th*  environment  chosen  is  a 
sever*  double  maxvelllan,  with  an  additional  component  representing  the  cold 
plasma  population.  The  results  presented  show  the  effects  of  small  changes  in 
th*  surfec*  characteristics  of  kapton  and  th*  incident  environment. 

2  THE  EQUIPOT  CODE 

The  development  of  EQUIPOT  is  based  on  s  number  of  requirements.  These  ate 
described  below,  together  with  details  of  the  code: 

(1)  The  baeio  geometry  ueed  by  the  oode  ehould  be  at  tittle  at  poetibla 
in  order  to  allow  modelling  of  differential  charging.  EQUIPOT  compute*  th# 
equilibrium  potential  of  s  patch  of  dielectric  material  (or  floating  conductor) 
on  an  undarlying  spherical,  conducting  spacecraft.  Tha  solar  illumination  of 
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tha  patch  and  the  structure  may  ba  varied  indapandantly .  The  Dabye  ahaath  which 
develop a  around  tha  patch  may  ba  aet  aa  either  planar  or  apherical.  The  simple 
geometry  uaad  by  EQUIPOT  la  shown  in  Fig  1.  It  ia  intended  that  the  planar,  or 
thin  ahaath  approximation  raaulting  from  ahaath  limited  collection  ia  used  in 
conjunction  with  LEO  environments,  where  tha  Debye  length  is  of  the  order  of  a 
few  millimetres;  the  spherical  or  thick  ahaath  approximation,  resulting  from 
orbit  limited  collection  ia  more  representative  of  the  conditions  which  prevail 
in  GEO  where  the  Debye  length  ia  of  tha  order  of  hundreds  of  metres,  much  larger 
than  a  typical  spacecraft.  Clearly,  these  two  approximations  represent  limiting 
cams,  realistic  conditions  will  fall  between  these  extremes,  aa  shown  in  Fig  2. 
Following  Ref  2,  the  currents  considered  are  due  to  incident  electrons,  back- 
scattered  electrons,  incident  ions,  secondary  electrons  from  incident  electrons 
and  ions,  photo-emission  and  conduction.  Tha  latter  involves  the  thickness  and 
conductivity  of  an  insulating  matarial,  or  cha  thicknasa  and  conductivity  of  a 
dielectric  substrate  for  an  isolatad  conductor.  Tha  ralative  permittivity  of  the 
dielectric  in  each  case  affects  the  capacitance  of  the  surface  element  and  the 
resulting  charging  time. 

(ii)  The  numerical  methods  used  throughout  the  code  should  be  as  simple 
as  possible  and  run  times  should  be  consistent  with  an  interactive  code.  EQUIPOT 
uses  a  simple  voltage  stepping  algorithm  in  crder  to  determine  the  equilibrium 
potential.  Mathematically ,  finding  the  equilibrium  potential  ia  a  straightforward 
root-finding  axarciae;  tharefore  it  would  be  possible  to  employ  fast  and  efficient 
algorithms  to  determine  Che  solution.  However,  simply  stepping  the  surface 
potential  from  an  arbitrary  starting  point  (usually  0  volts),  computing  tha  net 
current  at  aach  point  until  the  solution  ia  reaches  ensures  that  multiple  root 
cases  are  handled  properly,  and  also  shows  how  each  component  of  tha  net  current 
changes  aa  a  function  of  surface  potential.  EQUIPOT  also  includes  a  facility  to 
compute  the  current-voltage  -alationship  for  a  surface  between  user  defined 
potentials.  Fig  3  shows  Cha  aimpla  voltage  stepping  algorithm  used  by  EQUIPOT. 

(iii)  The  environment  definition  should  be  as  flexible  as  possible. 

EQUIPOT  accepts  plasma  environments  as  cables  of  measured  particle  spectra, 
single,  double  or  triple  maxvellian  components,  or  combinations  of  all  theta. 

For  example,  if  measured  spectra  are  available  for  a  given  energy  range,  ic  is 
poaaibla  to  append  maxvellians  which  giva  realistic  fluxes  outside  this  rsnga 
Altamativaly,  tha  environment  may  ba  defined  antirely  by  msxvellian  components. 
EQUIPOT  splits  tha  incident  parcicla  spectra  into  a  aeries  of  monoanergatic 
'swarms'  and  computes  ths  net  currant  due  to  aach  'swarm'.  Hence  tha  spasd  and 
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accuracy  of  the  simulation  ara  function*  of  tha  number  of  equal  logarithmic 
energy  ttep*  uaad  to  define  tha  incident  particle  energy  range  (usually  from 
1  eV  to  twenty  times  the  temperature  of  tha  hottest  maxwellian  component).  In 
general,  it  has  been  found  that  100  steps  per  integration  give  sufficiently 
accurate  and  quick  results. 

(iv)  Material  properties  should  be  defined  in  the  standard  NASCAP  form. 

•x 

Although  EQUIPOT  uses  the  same  general  format  as  NASCAP  ■  there  are  some 
important  differences.  The  properties  of  thickness,  conductivity  and  dielectric 
constant  ralata  to  the  insulating  patch,  or  to  the  dielectric  separating  the 
floating  conductor  from  the  spacecraft  structure.  Hence,  these  values  are  input 
to  the  cod*  separately  from  the  normal  material  properties.  Secondly,  NASCAP 
replaces  some  of  the  input  property  values  relating  to  secondary  electron  emission 
with  re-computed  values.  Since  EQUIPOT  treats  secondary  eleccron  emission  (SEE) 
differently,  this  method  is  not  adopted, 

(v)  Emphasis  should  be  placed  on  modelling  the  secondary  electron  emission 
yield  function.  This  plays  a  crucial  role  in  the  computation  of  equilibrium 

4 

potential  .  Hence  EQUIPOT  supports  a  number  of  different  empirical  and  theor¬ 
etical  secondary  electron  emission  yield  models:  these  will  be  referred  to  as 
the  'Kat2'  model1,  the  'Whipple/Dionne'  mode  1  6  and  Che  'Stcmglass'  model6. 

EQUIPOT  also  supports  normal  or  isotropic  SEE  yield  functions. 

(vi)  The  code  should  be  an  interactive ,  'engineering'  tool.  EQUIPOT  has 
been  developed  with  a  menu-driven  user  interface,  pre-defined  material  and 
environment  definition  files  and  a  continuously  displayed  panel  indicating 
current  set-up  ( eg  floating  and  structure  material  type,  integration  step  size, 
etc).  For  a  given  set  of  input  parameters,  the  equilibrium  potential  is  com¬ 
puted  in  a  time  of  typically  1-7  minutes,  allowing  many  parmutations  and  com¬ 
binations  of  parameters  to  be  assessed  at  a  tingle  session. 

3  SENSITIVITY  ANALYSIS 
3. 1  Method 

The  analysis  was  performed  assuming  a  kaption  patch  on  a  spherical, 
aluminium  spacacraft.  The  structure  is  assumed  to  be  sunlit,  and  therefore 
floats  at  a  few  volt*  positive.  The  surface  properties  of  kapton  ara  given 
below7: 
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Thickness 
relative  permittivity 
conductivity 
atomic  number 
Max.  normal  SEE  yield 
Energy  at  max  normal  SEE  yield 
yield  of  1  KaV  proton* 
energy  at  max  proton  yield 
phocoelactron  yield 


23  uo 
3.0 

1 . OE— 1 5  mho/m 

5.0 

1.9 

0.20  KaV 
0.4  55 
140.0  KeV 
2.0B-J  Am'2. 


The  sensitivity  analysii  computet  the  change  in  equilibrium  potential  as  each 
parameter  on  thit  liat  it  varied.  Where  appropriate,  the  input  parameters  are 
varied  by  10X  in  each  direction.  The  environment  chosen  for  the  study  is  a 

g 

variation  of  the  worst-case  SCATHA  environment  ,  modified  to  give  an  equilibrium 
potential  of  -10. 0  KV  when  applied  to  the  definition  of  kapton  given  above.  This 
provides  a  convenient  datum  which  is  used  throughout  the  study,  since  the 
percentage  change  in  the  result  is  simply  the  amount  that  the  equilibrium 
differs  from  -10. 0  KV  divided  by  100.  The  double  maxwellian  environment  with 
an  additional  cold  plasma  component  is  defined  as  follows: 


Ne, 

■ 

0.9 

CO 

cm  \ 

T*» 

■ 

600.0  «V 

Ne2 

m 

2.1 

T*2 

■ 

26000.0  *V 

h.3 

m 

0.1 

CO  3p 

t.3 

m 

1.0  *V 

m 

1 .0 

cm  \ 

-3 

Ti, 

m 

350.0  eV 

Ni2 

m 

0.7 

cm  » 
-3 

Ti2 

m 

25000.0  eV 

Ni3 

m 

0. 1 

cm  > 

Ti3 

1.0  *V  . 

This  differs  from  the  SCATHA  worst  case  environment  iu  two  ways;  firstly  a  low 
density,  cold  plasma  component  has  been  added  and  secondly  the  electron  and  ion 
densities  of  the  'hot'  component  have  been  adjuctad  to  give  the  required 
equilibrium  potential. 

It  should  be  noted  that  the  equilibrium  potential  of  -10.0  KV  obtained 
with  the  data  shown  above  used  the  Katz  model  of  secondary  electron  emission  and 
assumed  isotropic  particle  incidence.  The  time  taken  for  the  kapton  to  reach 
this  equilibrium  potential  was  4.0  hours. 

3.2  Variation  of  parameters  affecting  the  conductivity  currant 

The  thickness  of  kapton  and  it*  conductivity  affect  the  conduction  current 
to  the  structure  and  hence  the  equilibrium  potential.  Fig  4  show*  how  the 
equilibrium  potential,  Vg  varies  as  the  thickness  of  kapton  is  increased  from 
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5  urn  to  200  MO-  At  thicknesses  above  SO  pm,  V£  ia  relatively  insensitive  to 
changes  in  thickness,  whilst  below  this  value,  changes  in  thickness  cause  major 
changes  in  Vg  ,  Fig  S  shows  how  Vg  varies  for  25  pm  kapton  as  bulk  conduc¬ 
tivity  ranges  from  1  .0E-14  to  1.0E-16  mho/m.  The  significant  change  resulting 
froo  making  kapton  slightly  more  conductive  'i  evident. 

3.3  Variation  in  the  photoelectric  yield  current 

It  is  instructive  to  illuminate  the  kapton  patch  at  oblique  angles  of 
incidence.  EQUIPOT  allows  the  cosine  of  the  sun  angle  to  be  input,  where  1.0 
represents  normal  solar  incidence.  Fig  6  shows  how  Vg  is  affected  by  the 
variation  of  illumination  angle  from  complete  shadow  (90.0*)  to  near  72*  when 
negative  charging  is  completely  prevented  by  the  current  of  photoelectrons. 

3.4  Variation  of  SEE  current  for  isotropic  particle  incidence 

An  isotropic  flux  tends  to  increase  the  value  of  the  maximum  secondary 
electron  emission  yield  and  moves  its  position  to  a  higher  energy’’.  Two  of  the 
secondary  electron  emission  yield  models  supported  by  EQUIPOT  (the  Katz  and 
Whipple /Dionne  expressions)  include  an  expression  which  corrects  for  isotropic 
incidence.  The  positions  of  Lhe  peaks  in  the  corrected  models  are  shown  in 
Table  1  and  the  full  SEE  curves  are  shown  in  Fig  7.  Although  the  peak  positions 
are  in  good  agreement,  the  main  difference  between  the  curves  is  apparent  at 
high  incident  energies  where  the  Whipple/Dionne  expression  gives  a  much  higher 
yield. 

Teble  1 

Parameters  defining  the 
isotropic  SEE  yield  function 


SEE  model 

“Wx 

Wv) 

Katz  as  al,  1977 
Whipple,  1981 

2.74 

2.77 

338.0 

338.0 

rr 
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Table  2 


Rasul ci  of  variations  in 


the  isotropic  Katz  expression 


Var  iat  ion 

dMAX 

WeV) 

Default 

338  .0 

“max*  10x 

338.0 

dKAX  "  101 

2.466 

338.0 

W  101 

2.740 

371.8 

emax  '  101 

2.740 

304.2 

Table  2  shows  the  results  of  varying  the  maximum  yield  and  its  position 
for  the  (isotropic)  Katz  expression  by  10X  in  both  directions.  Since  only  the 
position  of  the  peak  is  being  changed  by  small  amounts,  and  the  shape  of  the  high 
energy  tail  of  the  distribution  changes  very  little,  the  equilibrium  potential 
is  affected  only  slightly  (leas  than  22). 

Table  3  shows  the  results  of  a  similar  analysis  performed  using  the 
Whipple/Dionne  expression.  Two  points  are  evident:  firstly  the  absolute 
equilibrium  potential  is  only  *2.8  kV,  significantly  different  from  the  results 
obtained  with  the  Katz  model,  and  secondly,  V£  is  much  more  sensitive  to  small 
changes  in  the  peak  position.  Once  again,  it  is  the  high  incident  energy  tail  of 
the  distribution  which  gives  a  such  higher  yield  than  the  Katz  expression,  and 
is  much  more  sensitive  to  changes  in  the  peak  position. 


Table  3 

Results  of  variations  in  the 


isotropic  Wh 


tmm 


Variation  dj^  Kj^ieV)  VE^ 
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It  it  not  intended  that  that*  ratulti  oak*  any  ttatcmcnt  about  the 
applicability  of  variout  SEE  models,  but  they  provide  a  clear  illustration 
that  the  choice  of  SEE  model  is  critical.  In  particular,  they  demonstrate  the 
importance  of  correctly  modelling  the  high  incident  energy  tail  of  the  SEE  yield 
function. 


3.5  Variation  in  SEE  current  for  normal  particle  incidence 

It  is  instructive  to  repeat  the  calculations  of  the  previous  section  for 
SEE  yield  functions  based  on  normal  particle  incidence.  It  should  be  noted  that 
the  expression  for  backscattered  electrons  is  also  corrected  for  normal  incidence. 
For  this  case,  EQUIPOT  also  supports  the  commonly  used  Sternglass  expression  as 
well  as  the  normal  versions  of  the  Katz  and  Whipple/Dionne  expressions.  Fig  8 
is  a  plot  of  the  three  functions.  By  definition,  all  three  curves  share  the  same 
peak  position  although  it  is  once  again  evident  that  the  main  difference  occurs 
at  high  incidenc  energies.  Table  4  shows  tha  results  of  mall  (101)  movements 
of  the  position  of  the  peak  on  equilibrium  potentials  obtained  with  each  of  the 
three  expressions.  For  the  Katz  yield  function,  changing  to  normal  incidence 
has  reduced  the  total  seconda  -y  yield,  such  that  the  equilibrium  potential  is 
now  near  -IS  kV,  Both  the  Katz  and  Sternglass  expressions  are  relatively 
insensitive  to  moving  the  location  of  the  peak,  in  fact,  small  changes  in  the 
energy  of  the  peak  have  almost  no  effect.  Changing  the  maximum  yield  by  101 
induces  changes  of  less  than  1Z  in  the  equ  . librium  potential.  V£  computed  with 
the  Whipple/D ionne  expression  changes  by  almost  UX  for  101  changes  in  the  position 
of  the  peak. 


Table  4 


Results  of  varying  che  normal  SEE  function 


Variation 

dMAX 

Wv) 

SEE  function 

Katz 

VE 

Whipple 

VE 

Sternglass 

VE 

Default 

1.90 

200.0 

-15008 

-11913 

-1564  2 

Vx  +  101 

2.09 

200.0 

-14941 

-11483 

-15640 

W  -  ,0* 

1.71 

200.0 

-15075 

-12337 

-15644 

‘max  ♦  '« 

1 .90 

22C.O 

-15007 

-11765 

-15637 

emax  -  102 

1.90 

180.0 

-15009 

-12068 

-15647 
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At  observed  for  isotropic  incidence,  the  equilibrium  potential  is  ouch  more 
sensitive  to  changing  the  secondary  yield  model  rather  than  changiug  the  peak 
position  for  any  single  model. 

3.6  Changes  in  the  current  of  beckscattered  electrons 

The  model  used  for  the  beckscattered  electron  yield  is  that  described  by 
Whipple^.  The  yield  is  a  function  of  the  atomic  number  of  the  surface  material 
end  also  whether  particles  are  incident  normally  or  isotropically .  EQUIPOT  uses 
the  isotropic  correction  automatically  if  the  SEE  yield  function  is  chosen  as 
isotropic,  hence  affects  of  the  different  backscatter  coefficient  for  normal  or 
isotropic  incidence  are  included  in  the  analysis  of  the  secondary  electron  yield 
function.  Therefore  the  atomic  number  provides  the  only  parameter  which  might 
sensibly  be  varied  over  a  limited  range  in  order  to  investigate  changes  in  the 
energy  dependent  backscatter  coefficient.  The  results  are  shown  in  Table  5  and 
demonstrate  that  the  effect  of  changing  2  is  small  (less  than  ll  for  102 
changes  in  Z  )  . 


Table  5 
Result  of 

small  changes  in  2 


- -  • 

Variation 

Z 

Yv) 

Default 

5.0 

-10000 

Z  ♦  101 

5.5 

-9933 

Z  -  102 

4.5 

-10077 

3 . 7  Changes  in  the  secondary  electron  current  due  to  on  impact 

The  yield  function  used  is  that  described  by  Whipple^  and  is  defined  by  two 
parameters;  the  yield  for  incident  I  KeV  protons,  and  the  energy  at  the  maximum 
yield.  The  results  of  varying  these  parameters  by  102  are  shown  in  Table  6.  All 
resulting  changes  in  V£  are  less  than  52.  However,  it  is  interesting  to  note 
the  effect  of  neglecting  this  current  component;  the  equilibrium  potential 
approaches  “2 1  KV. 
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Tabli  6 


Effact  of  varying  SEE  due  to  ion  impact 


Variation 

Yield 
(1  KaV) 

emax 

(KeV) 

vv> 

Default 

0.455 

140.0 

-10000 

Yield  at  1  KeV  *  10X 

0.501 

140.0 

-9604 

Yield  at  1  KeV  -  10X 

0.410 

140.0 

-10483 

Max  yield  energy  *  10X 

0.455 

154.0 

-9965 

Max  yield  energy  -  10X 

0.455 

126.0 

-10049 

No  SEE  due  to  ions 

-20797 

3.8  Chant* a  in  the  definition  of  the  plasma  environment 

Finally,  ic  ia  important  to  aiaaaa  tha  affect*  of  making  (mall  change*  in 
the  density  and  teogxrature  of  tha  maxvellian  component*  used  throughout  this 
study.  Firstly,  the  effect  of  the  cold  plasma  added  to  the  sever*  SCATHA 
environment  must  be  investigated.  A  density  of  0.1  cm'3  ia  too  small  to  be 
measured,  the  effects  of  both  removing  it  completely  and  also  increasing  it  are 
summarised  in  lable  7.  This  is  a  vary  significant  result  and  shows  how  tha 
presence  or  absence  of  cold  plasma  at  GEO  can  significantly  affect  spacecraft 
charging. 


Table  7 

Effect  of  varying  the  thermal  ion  density 


Variat ion 

Ni 

(cm  3) 

Ti 

(*V) 

vE(v) 

Thermal  ions  absent 

0.0 

1.0 

-16056 

Default 

0.1 

1 .0 

-10000 

Increased  density 

0.2 

1 .0 

-7829 

Table  8  shows  the  results  of  varying  tha  density  and  taoparature  of  the  two 
main  electron  and  ion  plasma  components  by  101.  The  colder  (600  aV)  electron 
component  does  not  contribute  to  the  current  balance  when  the  spacecraft  potential 
exceeds  several  thousand  volts  negative,  hence,  as  observed,  it  has  no  effect  on 
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a  spacecraf t  equilibrium  potential  of  around  -10  KV.  Changing  the  temperature 
and  density  of  the  'hot'  electron  component  reaulta  in  change*  of  leas  than  81 
in  Vg  ,  vhilet  aimilar  changes  in  the  'hot'  ion  component  give  rise  to  changes 
in  Vg  of  lees  chan  IX.  It  la  intareating  to  note  that  an  increase  in  the 
temperature  of  the  'hot'  ion  component  makes  Vg  leas  negative,  whilst  increasing 
the  temperature  of  the  cold  component  haa  the  opposite  affect,  although  both 
changea  are  small  (lasa  than  IX). 


Table  8 

Effect  of  varying 

£laamajdan*^t2_£Ili 

temperature 


Variation 

Vg(V> 

Default 

-10000 

Ne,  ♦  10X 

-10651 

Ne‘  -  10X 

-9362 

Ni,  ♦  10X 

-9925 

Nij-  I0X 

-10082 

Te,  ♦  10X 

-10721 

Tej  -  10X 

-9282 

Ti,  ♦  10Z 

-9975 

Ti,  -  10X 

-10032 

Ne  ♦  10X 

-10000 

Ns  -  10X 

-10000 

Ni.  ♦  10X 

-9818 

Ni]  -  10X 

-10225 

Te  ♦  10X 

-10000 

Te]  -  10X 

-10000 

Ti  ♦  10X 

-10085 

Ti]  -  10X 

-9911 

4  con;l’J5ION8 

Although  the  scope  of  this  study  is  limited  to  a  single  material  in  a  singla 
type  of  environment,  it  haa  yielded  some  valuable  results.  Two  main  features  hava 
emerged!  the  Importance  of  correctly  modelling  the  high  incident  energy  part  of 
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the  secondary  electron  emission  yield  function  and  the  importance  of  the  role 
played  by  a  cold  plasma  population. 

The  teeparature  of  the  'hot'  alectron  population  (23  KeV)  chosen  for  this 
study  corresponds  to  a  position  on  tha  SEE  yield  curve  of  between  75  and  125  times 
the  energy  of  the  yield  peak.  It  is  at  these  energies,  where  the  yield  is  less 
than  0.2,  where  various  modals  disagraa.  Further  SEE  yield  measurements  at  high 
incident  energies  are  necessary  to  rssolve  this  problem.  It  may  be  more  sensible 
for  codes  such  as  EQUIPOT  to  use  direct  SEE  yield  data,  rather  than  rely  on  an 
espirieal  relationship. 

Since  thermal  ions  are  rarely  completely  absent  from  the  geostationary 
plasma  environment,  particularly  during  prolonged  geosmgnetically  quiet  periods, 
it  is  essential  to  include  their  affect  in  any  charging  simulation. 

EQUIPOT  has  proved  to  be  a  vary  flexible  tool  in  performing  this  type  of 
analysis.  Future  studies  are  proposed  which  include  investigation  of  the 
behaviour  of  materials  in  measured  geostationary  plasma  environments  with  a  view 
to  explanations  of  the  'eclipse'  and  'barrier'  charging  events  observed  on  the 

9 

Meteosat-2  spacecraft  . 
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Fig  4  Equilibrium  potential  as  a  function  of  kapton  thickness 


Tig  6  Equilibrium  potential  as  a  function  of  illumination  angle  for  a  kapton  surface 
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Fig  7  Secondary  electron  yield  for  kapton;  isotropic  particle  incidence 


fig  8  Secondary  electron  yield  for  kapton;  normal  particle  incidence 
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